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It is shown in the framework of the fully relativistic Dirac-Fock treatment of 
photoionization and radiative recombination processes that taking into account all 



, significant multipoles of the radiative field is of considerable importance at electron 



energy higher than several keV. For the first time, we show that the relativistic 
Maxwell-Bolzmann distribution of continuum electrons should be used in hot thermal 
plasmas. This decreases the radiative recombination rate coefficient up to several 
multiplies compared to the non-relativistic distribution commonly used. 
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■ The photoionization and radiative recombination cross sections as well as the radiative recom- 

o : 

bination rate coefficients are required for estimates of ionization equilibria and thermal balance 
in terrestrial and astrophysical plasmas contaminated by various ions. In fusion reactors at tem- 
^ ■ peratures above several keV, impurity atoms of various elements may be stripped to bare nuclei. 

At temperatures around 1000 keV, tungsten atoms are fully stripped In astrophysical ob- 
jects such as stellar black-hole binaries and Seyferrt galaxies, the plasma temperature may reach 



150 keV 



2]- 



At sufficiently high kinetic electron energy Ek and for highly charged ions, multipole and rela- 
tivistic effects should be taken into account in calculations of photoionization cross sections (PCS), 
radiative recombination cross sections (RRCS) and radiative recombination rate coefficients. The 
effects were considered beginning with pioneering works . Nevertheless, these effects are 
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usually neglected in the application of these processes in plasmas (see Refs. in 6|, |7|). The most 
extensive advanced calculations by Badnell [7] were performed using the electric dipole and semi- 
relativistic approximations for electron energies to 1.36Z 2 keV and for Z < 54, that is up to 
~ 4 MeV. 

We used the fully relativistic treatment of the photoionization process in recent calculations of 
PCS and RRCS for 31 ions of elements from the range 26 < Z < 74 [6]. Electron wave functions 
were generated in the self-consistent Dirac-Fock (DF) framework. We took into account all sig- 
nificant multipole orders of the radiative field. Previously, we performed relativistic calculations 
of total and differential RRCS for recombination of an electron with the H-, He- and Li-like ura- 
nium ions using this model with regard to the Breit electron interaction and the main quantum 



electrodynamic corrections 



The influence of the multipole effects was also considered in our 



studies of the total and differential PCS 



Exact relativistic benchmark calculations of RRCS including all significant multipoles were 
carried out by Ichihara and Eichler 12J for radiative recombination of the K, L and M electrons 
with bare nuclei with charge numbers 1 < Z < 112 . The results for a few representative cases 
were compared with those derived from the widely used non-relativistic dipole approximation in 
order to assess the accuracy of the latter. 

In the present paper, we discuss the influence of multipole effects on PCS, RRCS and rate 
coefficients as well as the influence of the relativistic effects such as the relativistic transformation 
coefficient between PCS and RRCS and the relativistic correction factor for the non-relativistic 
rate coefficient. This temperature depend factor is reported for the first time in this paper. 

The relativistic PCS in the i-th subshell per one electron can be written in the form 

Air 2 a 



W 
ph 



V£ (2L + 1)Q 2 ll (k) + LQl +1L (n) 
+ (L + - 2JL(L + 1) Q L -il{k)Q l+1l {k) 



(1) 



Here k is the photon energy in tuqC 2 , L is the multipolarity of the radiative field, k = (£—j)(2j+l), 
£ and j are the orbital and total angular momenta of the electron, a is the fine structure constant 
and Qal{ k ) is the reduced matrix element (for detailed expressions see jfj]). 

The cross section of the recombination process with the capture of an electron with energy 
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to the i-th subshell of the ion is expressed in terms of the corresponding PCS as follows 

a« = A qi c$ , (2) 

where qi is the number of vacancies in the i-th subshell prior to recombination. The transformation 
coefficient A can be derived from the principle of the detailed balance. The exact relativistic 
expression for the coefficient is written as Q, Q 

. k 2 ~ E k 

A ve \ = ~ ~ , E k = -. (3) 

2E k + El m c 2 

However in the majority of the RRCS calculations, the coefficient is used in the form 

k 2 

Anrel = - jr— , k = m C 2 k , (4) 

2m c 2 E k 

which may be obtained in the non- relativistic approximation from Eq. ([3]). The difference between 
<r„ obtained with Eq. ([3]) and Eq. (@J depends only on the electron kinetic energy E k and can be 
written as 



A nr cl — A vc \ ^ E k . . 

A iel ~ 2m c 2 

with the difference ~ 5% at E k = 50 keV and reaching ~ 100% at E k = 1000 keV. Consequently, 
at high electron energy, the relativistic expression ([3]) should be used in the RRCS calculations. 

The relativistic recombination rate coefficients a^(T) can be calculated using the thermal 
average over RRCS. In the present paper, the continuum electrons are described by the relativistic 



Maxwell-Boltzmann distribution function f(E) normalized to unity as follows 13] 



Here E is the total electron energy in units of ttiqC 2 including the rest energy, 9 = kpT/moC 2 is the 
characteristic dimensionless temperature, T is the temperature and kp is the Bolzmann constant. 
The function K2 denotes the modified Bessel function of the second order. The relativistic rate 
coefficient may be written as 

o£J(T) =< va® >= F iel (6) ■ a«(T), (7) 
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where v = (p/E)c is the electron velocity with the momentum p = \fE 2 — 1 and ay'iT) is the 
usual non-relativistic rate coefficient [14] 

oo 

« W (T) = (2/7r) 1 / 2 c- 2 (m k^T)- 3 / 2 q l J k 2 a^(k)e^- k ^ T Uk, (8) 

Si 

where is the photon energy and £j is the binding energy of the i-th shell. In Eq. (7), F re \(9) is 
the relativistic factor 

Frei(e) = \J\0 / K 2 {l/6)e^ e . (9) 

I — I 

Using the asymptotic expansion of the Bessel function K 2 (l/9) at large 1/9 [l5|], that is at low 
temperature, we arrive at the factor F Te \(9), an approximation to F Te \(9) 

*-w = 1 /( 1 + T , + ii J> + -)- (10) 

Eq. (JTOl) provides an excellent approximation for F TC \(9) with the terms through order 9 2 at 9 < 1. 
The factors -Frei^) and F re i(^) are compared in Fig. 1. The solid curve refers to the exact factor 
(Eq. (jHl)) and the dashed curve refers to the approximate factor (Eq. (TTD]) ). As can be seen, there 
is little difference between the two curves, the relative error is ~4% at kpT = 500 keV and 25% 
at k a T = 1000 keV. 
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FIG. 1: Exact factor F Te \(T) (solid) and approximate factor F TC \(T) (dashed). 

As is seen from Fig. 1, the use of the relativistic distribution instead of non-relativistic one 
results in a decrease of rate coefficients by a factor of 1.2 at plasma temperature kpT = 50 keV 
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and up to factor of 7 at kpT = 1 MeV. 

Let us next consider the influence of the multipole effects. The electric dipole approximation 
takes into account only terms with L—l in Eq. (1). As is well known, the dipole approximation 
holds at a low electron energy E k but breaks down at a higher energy. 
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FIG. 2: Subshell RRCS (in barns) calculated taking into account all multipoles L (solid) and in the 
dipole approximation (dashed). 



In Fig. 2, we compare RRCS obtained in the dipole approximation cr rr (dip) (dashed curves) 
with RRCS calculated with all multipoles <7 rr (L) making significant contribution (solid curves) for 
bare nuclei of two representative elements Fe (Z—26) and W (Z—74). The energy range under 
consideration is 1 keV < E k < 1000 keV. As is seen in Fig. 2, the curves begin to diverge noticeably 
even at several keV. At the highest energy 1000 keV, in the case of W 74+ , <r rr (dip) is smaller than 
the exact value o" rr (L) by a factor of ~5 for the Is shell and by a factor of ~40 for the 4/ 5 / 2 
subshell. Our calculations showed that the relative difference between the exact calculation of 
RRCS and the dipole approximation 

a rr (L) - cr rr (dip) 



mRCS 



(T rr (L) 



100% 



;n) 
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varies in the range ~3-20% for shells with different orbital momenta ii at E^—10 keV, ~ 15-50% 
at E k =50 keV and reaches several multiples at ^=1000 keV. The dependence of A RRCS on £j is 
shown to be considerable, A RR cs being larger with increasing £ f . The difference A RRCS was found 
to increase with Z, especially for high energy. 

Calculations show that to achieve accuracy ~0.01% in the PCS, e.g., for the Is and 3c? shells 
of the ion W 73+ one has to take into account all terms in Eq. (1) up to L=5 and L—8, respectively 
at E k =10 keV, to L=7 and 11 at E k =W0 keV and to L=19 and 31 at £ fc =1000 keV. 



Table I. Comparison of our PCS with results by Badnell [7J for the Is shell of the H-like ion 



Xe 53+ . A 



PCS 



Wph (present) — <7 p h( Badnell) }/cr ph (present) 



• 100%. 



<Tp h , Mb 



E k , keV Badnell Present Apes, % 
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1.318( 


-7) 
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3967. 
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350 (• 


-9) 


4.117( 
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In Table I, we compare our present PCS calculations with the corresponding results of Bad- 
nell j3] for the Is shell of the H-like ion Xe 53+ . The case of the one-electron ion is particularly 
convenient for checking the influence of the higher multipoles and the method of calculation be- 
cause there are no any inter-electron interactions. In this case, the PCS must be independent of 
the gauge used in calculations for correct wave functions. 

We found that our calculation is in excellent agreement with values from [12] where all mul- 
tipoles L were involved. For the Is shell of the H-like ion Xe 53+ , the two calculations coincide 
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with an accuracy of the three significant digits presented in [12J in the wide energy range 1 eV 
< Ek < 6000 keV. By contrast, PCS obtained by Badnell exceed our values and results from 12| 
by ~ 16% in the energy range < 4 keV and diminish progressively at higher energies becoming 
lower by a factor of ~ 8 at Ek ~ 1800 keV and a factor of ~ 30 at Eu « 4000 keV compared with 
our values. The comparison of our PSC values and results from 12j with calculation by Badnell 
[7] for the lighter ion Fe 23+ reveals a similar tendency, but smaller in magnitude. 

The reason of the difference at low energies is unclear for us because the non-dipole terms 
make a small contribution at low energies (see Fig. 2). It is possible that the difference arises from 
the methods of calculation used in [j}. The difference at high energies (> 100 keV) must be due 
to neglect of the higher multipoles and possibly also due to the semi-relativistic approximation 
adopted in 



30 
20 

10 


40 
g 30 
<? 20 

10 



40 
30 
20 
10 


10 7 1 8 1 9 1 10 

Temperature (K) 

FIG. 3: Difference A a between rate coefficients calculated with using all multipoles and in the dipole 
approximation for the 2s (solid), 2p 1 / 2 (dashed), and 3d 3 / 2 (dash-dotted) shells. 

From the discussion above, it would be expected that the dipole approximation would also fail 
in calculations of rate coefficients at a high temperature T. In Fig. 3, we present the difference A Q 
between the exact a^(L) and the dipole a^(dip) values of partial rate coefficients. The difference 
is defined in the same way as in Eq. ( fTTj) . The difference A a is given for the 2s, 2py% and 3ds/2 
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electrons recombining with the He-like ions Fe 24+ , Xe 52+ and W 72+ . These shells are the lowest 
ones making a large contribution to the total rate coefficients. As is evident from the figure, the 
difference A Q is larger for the heavy, highly charged ions. The inclusion of higher multipoles may 
change partial rate coefficients by ~ 7% at temperature T = 10 8 K, by ~ 20% at T = 10 9 K and by 
~50% at T = 10 10 K for W 72+ . This means that total rate coefficients obtained within the dipole 
approximation have to be considerably smaller than the accurate values obtained with regard to 
all multipoles L. 

In conclusion we have clearly demonstrated the importance of multipole effects to the PCS, 
RRCS and rate coefficient calculations at electron energy of the order of 10 keV and higher. We 
have also showed that in hot plasmas, the relativistic Maxwell-Bolzmann distribution of continuum 
electrons must be used in the rate coefficient calculations. It should be noted that the relativistic 
rates may be similarly obtained for processes of dielectronic recombination and the electron impact 
ionization in hot plasmas. 
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